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"on ! Abstract 
O 

We study the dynamics of a probe soliton propagating in an optical fiber and exchanging energy in 

CN 

fast collisions with a random sequence of pump solitons. The energy exchange is induced by Raman 

2: 

scattering or by cubic nonlinear loss/gain. We show that the equation describing the dynamics 
of the probe soliton's amplitude has the same form as the equation for the local space average of 
energy dissipation in random cascade models in turbulence. We characterize the statistics of the 
probe soliton's amplitude by the r q exponents from multifractal theory and by the Cramer function 
S{x). We find that the nth moment of the two-time correlation function and the bit-error-rate 
contribution from amplitude decay exhibit power-law behavior as functions of propagation distance, 
^ ■ where the exponents can be expressed in terms of r q or S(x). 
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The dynamics of localized patterns in the presence of noise and nonlinear effects is a rich 
subject that is of major importance in many fields, including solid state physics jlj, turbu- 
lence and optics Fiber optics communication systems, which employ optical pulses 
to represent information bits, provide an excellent example for systems where noise and non- 
linearities have an important impact on pattern dynamics (3|. It is well established that the 
parameters characterizin g op tical pulses in fiber optics networks can exhibit non-Gaussian 



statistics 



10t ] . Yet, since optical fiber systems are only weakly nonlinear, it 



was commonly believed that the statistics of optical pulse parameters is very different from 
the statistics encountered in strongly nonlinear systems, such as turbulence and chaotic flow, 
where intermittent dynamics exists. Two recent studies of pulse propagation in multichan- 
nel optical fiber transmission systems in the presence of delayed Raman response obtained 
results that stand in sharp contrast to this common belief ll|, ll2j. Taking into account the 



interplay between Raman-induced energy exchange in pulse collisions and bit pattern ran- 
domness it was shown that the pulse parameters exhibit intermittent dynamics in the sense 
that their normalized moments grow exponentially with propagation distance. Furthermore, 
it was found that this intermittent dynamic behavior has important practical consequences, 
by leading to relatively large values of the bit-error-rate (BER), which is the probability for 
an error at the output of the fi 



The studies in Refs. 



12 



iber line. 

raise many intriguing questions, which are of fundamental 
importance in nonlinear optics, statistical physics, and chaos and turbulence theories. The 
main questions are: (1) Is the similarity between the fiber optics system and the turbulent 
one coincidental, or is it a consequence of a common underlying mechanism? (2) Is Raman 
scattering the only nonlinear process leading to intermittent dynamics of optical solitons? 
(3) Can the statistics of pulse parameters in the weakly nonlinear optical fiber system be 
analyzed by statistical mechanics tools that are used in the analysis of strongly nonlinear 
systems? (4) If the answer to (3) is yes, what predictions can be made for the main ob- 
servables using these statistical mechanics tools? In this Letter we address these questions 
in detail. We focus attention on propagation of a probe soliton undergoing many collisions 
with a random sequence of pump solitons in the presence of delayed Raman response. First, 
we show that for certain setups of the pump soliton sequence, the dynamic equation for 
the probe soliton's amplitude has the same form as the equation for the local space aver- 
age of energy dissipation in random cascade models in turbulence. The latter models play 



2 



a pivotal role in the analysis of large fluctuations of the local energy dissipation and in 
obtaining corrections to the scaling laws of the classical Kolmogorov theory of turbulence 



B, Q, Q, Q, Q, 0, Q. 
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201 ]. Thus, our study provides a surprising answer to question 



(1): the similarity between the weakly nonlinear system and the strongly nonlinear one is 
not coincidental. Second, we demonstrate that Raman scattering is not the only nonlin- 
ear effect leading to intermittent dynamics and that similar behavior should be observed 
in systems where the collision-induced energy exchange is due to cubic nonlinear loss/gain. 
In this sense our results are quite general, since they hold for the three main first order 
perturbations to the nonlinear Schrodinger (NLS) equation that lead to collision-induced 
energy exchange. Third, we obtain a positive answer to question (3) by showing that the 
statistics of the probe soliton's amplitude can be described by the r q exponents, which are 
commonly employed for analyzing multifractals and strange attractors in chaotic systems 
jltl 17, 18], and by relating the r q exponents to the Cramer function S(x) via a Legendre 



transform. Fourth, we examine the implications of this dynamic behavior on two major 
observables of the probe soliton: the nth moment of the two-time equal-distance correlation 
function and the BER contribution from amplitude decay. We find that both observables 
exhibit power-law behavior as functions of propagation distance, where the exponents can 
be expressed in terms of r q or S(x). We also reveal an intriguing similarity between BER 
dynamics and the dynamics of the <i-measure in coarsening of geometrical multifractals. 
It is interesting to note that several previous works studied the emergence of optical 



turbulence in nonlinear cavity ring resonators 21 
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241 ] . In one of these works (Ref. 



24j) the emergence of turbulent-like behavior was associated with the excitation of many 
strongly interacting spatio-temporal structures, but the phenomenon was not analyzed quan- 
titatively. In addition, direct mathematical links between dynamics in systems described 
by NLS and related models and dynamics observed in models of turbulent flow have been 



obtained in Refs. 



25 



26 



27|, 



281 ] (see also references therein). However, the propagation 



equations in these studies did not take into account the effects of Raman scattering, and the 
setups considered were different from the setups considered here. Furthermore, the statistics 
of the physical observables was not analyzed in terms of the random cascade model and the 
multifractal formalism that are employed in the current paper. 

Propagation of pulses of light in an optical fiber in the presence of delayed Raman response 



3 



is described by the following perturbed NLS equation 

id z V + d 2 t ^ + 2|^| 2 ^ = -e R tyd t \ty\ 2 , 



(1) 



where \I/ is the envelope of the electric field, z is propagation distance and t is retarded 
time. The term — e#\l/<9 t |\l/P accounts for the effects of delayed Raman response and e R 
is the Raman coefficient 29j. When e# = 0, the solution of Eq. ([T]) corresponding to 
a single soliton with frequency (3 is described by ^p{t : z) = rj^ exp(ixp) cosh -1 ^), where 
xp — VP {t ~ Up ~ 2/3 z), xp — a p + Pit ~ Up) + {Vp ~ P 2 ) z i an d f]pi cap and y@ are the soliton 
amplitude, phase and position. 

Consider a single collision between a probe soliton with frequency ft — and a pump 
soliton with frequency /3. In a fast collision |/3| ^> 1. Assuming in addition that e# <C 1, 
the main effect of delayed Raman response on the collision is an 0(€r) change in the soliton 



amplitude 



1Q,& 



31 



32|: 



(2) 



10, B 



Ar] = 2r] r] f3 sgn(P)e R . 

The effect of the collision in order €r/ /3 is a frequency shift given by 
— (8?7o?7/3ei?)/(3|/3|). This effect and effects of order e 2 R and higher can be neglected for the 
dynamical setups considered below. In addition, for other types of perturbations, such as 
those due to third order dispersion, the collision-induced changes in amplitude and frequency 
are of higher order in both the parameter e characterizing the perturbative process and l/\/3\ 



(see, e.g., Refs. 
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34 
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We now describe propagation of a probe soliton under many collisions with a random 
sequence of pump solitons. The pump solitons are located at time slot centers, and each time 
slot can be either occupied or empty. The occupation state of the jth time slot is described 
by the random variable Q: Q = 1 with probability s if the slot is occupied and Q = 
with probability 1 — s otherwise. It is assumed that different time slots are uncorrelated: 
(dCj) = s 2 if i 7^ j. The probe soliton is initially located at yo(0) = and the frequencies of 
the pump solitons are f3j = A/3 > 0. Since we look for power-law behavior of the physical 
observables we assume that the pump solitons are initially located at ypj(0) = — a J T, where 
T and a > 1 are constants. Therefore, the collision between the jth pump soliton and the 
probe soliton occurs at a distance Zj, given by: Zj = a^Azc, where Az^ = T/(2A/3). 

Using Eq. (T2]) and assuming that %j(0) = ^(O) for all pump solitons we find that the 
probe soliton's amplitude after J collisions is r) (zj) = rj (zj-i)[l + 2e^?7^(0)O]; leading 



to t] (zj) = Vo(0)Tlj=iWj, where Wj = 1 + 2€Rrjp(0)Q. Compensation of average cross 
talk effects can be introduced in a straightforward manner. In this case the probe soliton's 
amplitude after J collisions is rj (zj) = r] (zj_i)[l + 2en7}p(Q)(Q — s)], resulting in 

j 

Vo(zj)=Vo(0)l[W 3 , (3) 

j=i 

where Wj = l + 2ejir}p(0)(Q — s). Notice that (Wj) = 1, and therefore (rj (zj)) = 1. Equation 
([3]) has the same form as the equation for the local space average of energy dissipation in 
random cascade models in turbulence. (Compare with Ref. Q], p. 166). In this equivalence, 
Zj, Azi 1 ^ and r] (zj) } play the roles of eddy size I, upper turbulence cutoff l Q , and energy 
dissipation of eddies of size I, ej, respectively. The evolution of both rj and ei is multiplicative 
and dissipative. Notice that in the fiber optics system energy cascades between successive 
collisions, whereas in the turbulent model the cascade is from large eddies to smaller ones. 
Using Eq. (j3J) and the statistical independence of the Wj factors, we obtain 

Wzj)) = « (zj/AzMf ^ . (4) 

Q, 



We define the r q moments in a similar manner to Refs. 



13, 



18|]: r q = log a (W q ). For the 



fiber optics system: 

(W q ) = sw\ + (1 - s)w%, (5) 

where w\ = 1 + 2(1 — 5)6^77/3(0) and W2 = 1 — 2se/j77/3(0). The r q curve obtained by using 
Eq. © is plotted in Fig. fflfor e fi = 0.03 (f = 0.2 ps), a = 1.25, ^(0) = 1, and different 
s- values. 

We emphasize that similar dynamics of the probe soliton's amplitude is expected in 
systems described by perturbed NLS equations, where the perturbation is due to cubic 
nonlinear loss/gain. For these systems the — e^dt \^\ 2 term on the right hand side of Eq. 
([I]) is replaced by =Fe c | \I/ 1 2 \I/, respectively, where e c is the cubic nonlinear loss/gain. The main 
effect of a fast collision in this case is an amplitude change, which is given by an equation of 
the form (j2J) with sgn(/3)e# replaced by =F2e c /|/3|. Thus, the results of this Letter are quite 
general since they hold for Raman scattering and nonlinear loss/gain, which are the three 
main first order perturbations to the NLS equation that lead to energy exchange in pulse 
collisions. 
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FIG. 1: The r q curve for the probe soliton's amplitude for er = 0.03 (pulse width = 0.2 ps), 
a = 1.25, and 77/3(0) = 1. The solid, dashed and dotted lines represent the values obtained in the 
case where 77/37(0) is deterministic with s = 0.5, s = 0.1, and s = 0.9, respectively. The dashed- 
dotted, dashed-dotted-dotted, and short-dashed curves correspond to r q in the case where 77/37(0) 
is random with s = 0.5, Arjpj = 0.5, and p = 0.5, p = 0.1, p = 0.9, respectively. 



It is possible to relate the r q exponents to the Cramer function S(x), characterizing the 
statistics of the pump soliton bit pattern, by a method similar to the one used for random 
cascade models and multifractal sets [ly, Q- For this purpose we notice that the 
variable m = ^2j=x Cj * s binomially distributed with probability density function (PDF) 
P(m; J) = [J!s m (l — s) J ~ m ]/[m\(J — m)!]. Using Stirling's formula we obtain 

P(m; J) ~ [S"(x)} 1/2 exp[JS(x)]/ (2ttJ) 1/2 , (6) 

where x = m/J and 

S(x) = -x ln(x/s) - (1 - x) ln[(l - x)/(l - s)]. (7) 

To obtain the relation between r q and S(x) we notice that a jTq = (Ylj=iWj), and conse- 
quently: a jTq = (exp ^qJ2j=i m [1 + ^ e R(Cj ~ Expressing the sum in the exponent's 
argument in terms of x, w± and W2 and using the large-deviations theorem we obtain 

a Jr q ~ /da; — exp {qj [xln(wi/w2) 

J L 2?r J 

+ ln(w 2 )] + JS(x)}. (8) 

Employing a saddle point approximation we obtain the Legendre transform relating r q and 
S(x): 

T q = {q [x\n(w l /w 2 ) + ln(io 2 )] + S(x)} / lna, (9) 

6 



where x(q) is determined by the conditions gln(wi/w 2 ) = — S'(x), and S"(x) < 0. 

Since e«<Cl the Wj factors in the fiber optics system are only slightly different from 1. 
It is therefore interesting to look for ways to enhance the randomness of the collision-induced 
energy exchange. Here we describe one disorder enhancement mechanism, which also leads 
to a natural generalization of the model described by Eq. (J3]). Consider a setup in which the 
amplitude of the jth pump soliton is given by 77/3? (0) = 77/3(0) ± Arjp with probabilities p and 
1 — p, respectively. A straightforward analysis shows that the 770 dynamics is still described 
by Eq. ([3]). However, the Wj factors can now attain three values: Wi = l + 2eij[(l — 5)77^(0) + 
(1 — 2p)sArj[3 + Arip] with probability ps, W2 = 1 + 2cr[(1 — s^^O) + (1 — 2p)sArjp — Ar)p] 
with probability (1 — p)s, and w 3 = 1 + 2eR[— 37^(0) + (1 — 2p)sAr] /3 ] with probability 1 — s. 
The average of W q is 

(W q ) = psw q 1 + (l-p)sw q 2 + (l-s)wl (10) 

The T q curve for this model is shown in Fig. [U One can see that r q attains larger values when 
r]/3j(0) is random for s = 0.5, p = 0.5, and p = 0.9, compared with the case where 77^(0) is 
deterministic and s = 0.5. Thus, the cumulative effect of energy exchange in the collisions 
is indeed enhanced by randomness of ?7/3j(0). To obtain the Legendre transform between r q 
and S(x), we denote by m^i = 1,2,3, the number of occurrences of pump solitons with 
Wj = Wi. The PDF of mi and 7712 is trinomial: Pimx^m^ J) = {J\(ps) mi [(l — p)s] m2 (l — 
s) J ~ mi ~ m2 }/[mi\m2\(J — mi — 7712)!]- Using Stirling's formula we arrive at 

u / n J 1/2 exp[JS(x u x 2 )] 

P( mi ,m 2 ; J) ~ — -— , (11) 

2tt [m 1 m 2 (J — rrii — m 2 )\ 

where x\p = mi t2 /J, and 

S(xi,x 2 ) = -x 1 \n[x 1 /(ps)] - a; 2 ln{x 2 /[(l - p)s)} 

_(1 _ Xl - X2 ) l n [(i _ Xl - X2 )/(i _ s )]. (12) 

By employing the large-deviations theorem and a saddle point calculation we obtain the 
following generalization of the Legendre transform given by Eq. Q: 

r q = {q [xi ln(wi/K; 3 ) + x 2 ln(w 2 /w 3 ) + ln(tt> 3 )] 

+S(x 1 , x 2 )\ /In a, (13) 
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where xi(q) and X2(q) are determined by 

OS 

q\n(w i /w 3 ) = -—\ Xi=xM} 2 = 1,2, (14) 

d 2 S/d 2 xi < 0, and (d 2 S / d 2 x\)(d 2 S / d 2 x 2 ) — (d 2 S/dxidx 2 ) 2 > 0. Since the exact expression 
for (W q ) [Eq. (flQl) ] can be retrieved by using Eqs. (fT2l - ([Tll) we conclude that the r q 
exponents give a correct characterization of the statistics of rj . In addition, Eqs. (fT2l - f[T4"j) 
can be further generalized for the case where the Wj factors attain any finite number of 
values. 

We now discuss implications of the dynamics of the probe soliton's amplitude. We start 
by considering the nth moment of the two-time equal-distance correlation: 

C& = ((^oor/oi) n )/((^oo)"(r/oi) n ). (15) 

The function Cq" measures correlation between the amplitudes 7700 and 7701 of two probe 
solitons, whose initial positions are ?/oo(0) = and yoi(0) = AT, respectively, where A is 
a constant. We remark that high moments of velocity correlation function along a given 
direction play an important role in turbulence theory 2J. From Eq. ([3]) it follows that 
([Voo(zj)voi(zj)) n ) = ^oo(0)^oi(0)(ni=i%nj=r%i), where Jmax is the number of col- 
lisions experienced by the 01 probe soliton, and Wji = Wj for j = 1, . . . , j max . Using 
this relation, the statistical independence of the Wjo factors and the definition of the r q 
exponents, we obtain 

Tin |~ "1 T n 

(16) 



Zj 



zj - Z X \ 



where Z\\ = AAz^\ Thus, the nth moment of the two-time correlation function is a product 
of power-laws of two different scaled distances. 

One of the most important quantities characterizing the performance of fiber optics sys- 
tems is the BER. Since the BER is often determined by the tail of the PDF of the pulse 
parameters one can expect that it would be closely related to the Cramer function. We show 
that this is indeed the case for the system described here. We consider the setup where the 
amplitudes of the pump solitons are deterministic [^j(O) = ^(O)] and focus attention on 
the contribution to the probe soliton's BER due to amplitude decay, BER^. This contribu- 
tion is defined by: BER,, = JJ 7 " 1 dr] F(r]o), where F(r} ) is the amplitude PDF and r\ t h is the 
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threshold for an error. Using r]o(zj) = t/q^w^w^ m and Eq. (J6]) we obtain 

BER 



g jWexplJSim/J)] 



(17) 



m=0 [2vrm(J-m)] 1/2 

where m t h is the solution of rjo(m) = 77^. Since S(x) is an increasing function of x the main 
contribution to the sum on the right hand side of Eq. (TTTI) comes from a close neighborhood 
of ra t h- Taking into account the two leading terms we arrive at: 



BER™ 





~-S"{x)~ 


1/2 




Si x th) 
In a 


1 + a lna 




2nJ 









(18) 



where Xth = mth/ J- Thus, in the leading order BER^ grows like a power law with propaga- 
tion distance, where the exponent is the value of the Cramer function at the error threshold. 
The eR-dependence of BER,, in the range 0.0075 < e R < 0.012 (0.5ps < f < 0.8ps) is shown 
in Fig. [2] for 770(0) = 1, 7^(0) = 2, r, th = 0.7 and J = 27. For the choice a = 1.25, T = 5, 
A/5 = 40, and 02 = — 4ps 2 /km, for example, we obtain 7.96THz < Au < 12.74THz for the 
frequency difference, and 1.07km < X27 < 2.75km for the propagation distance. It is seen 
that both Eq. (TlBl) and the first term on the right hand side of Eq. (fl~8l) are good approx- 
imations to the exact result over a wide range of BER^ values. It is interesting that the 
dynamic behavior of BER^ is very similar to the behavior of the <i-measure of geometrical 
multifractals during the late stage of coarsening (compare Eq. f[T7|) with Eq. (9) in Ref. 



36|) 



In many cases the dynamics of the other three soliton parameters is coupled to the am- 
plitude dynamics and as a result, these parameters can be strongly influenced by amplitude 
fluctuations. For the system considered here the most important effect is due to the Raman- 



induced self frequency shift, which is given by 
frequency shift leads to a position shift: 



231 : fo{z) = -(8e R /15) J Q Z dz'^(z'). This 



Vo{z) 



lQe R 
15 



dz' / dA V) 
u Jo 



(19) 



that can give significant contribution to the probe soliton's total BER. It is therefore im- 
portant to evaluate the impact of this process on the total BER in comparison with the 
contribution BER^ coming solely from amplitude decay. Since we do not have an analytic 
expression for the PDF of y we carry out Monte Carlo simulations with Eqs. ([3]) and ([191 . 
We define the relative position shift yo = Ho — (no) , where the average (7/0) is assumed to be 
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FIG. 2: Total BER and BER contribution due to pulse decay, BER^, vs the Raman coefficient €r 
for ?7o(0) = 1, rjp(0) = 2, a = 1.25, T = 5, and 02 = — 4ps 2 /km. The squares, up triangles, and 
circles stand for the exact result for BER^, the approximate expression given by Eq. (|18p . and 
the first term on the right hand side of Eq. (|18p . respectively. The down triangles and diamonds 
represent the total BER for A/3 = 40 and A/3 = 20, respectively. 

compensated by niters. For each realization of the Wj factors we compute the total energy 
at the detector at distance zy. 

I{zj) = r) {zj) {tanh[r /0 (z J )(T/2 - y ( Zj ))) + taii%o(>j)(T/2 + y {zj))} . (20) 

An occupied time slot is considered to be in error, if I(zj) < I t h = ^Vth tanh(^T/2). We 
use the same parameter values as described in the previous paragraph, but with two different 
values of A/3: A/3 = 40 and A/3 = 20 {38]. The total BER obtained in the simulations is 
shown in Fig. [2j It is seen that for A/3 = 40 the total BER is very close to BER^, that is, 
error generation is dominated by amplitude decay. In contrast, for A/3 = 20, error generation 
is dominated by the Raman-induced position shift, and as a result the total BER is much 
larger than BER,^. These results can be explained by noting that the smaller inter-collision 
distances for A/3 = 40 lead to a smaller total propagation distance, and consequently, the 
position shift is relatively small. For A/3 = 20, the inter-collision distances and the total 
propagation distance are large, resulting in relatively large position shifts. 

In summary, we studied the amplitude dynamics of a probe NLS soliton, exchanging 
energy in fast collisions with a random sequence of pump solitons. We showed that the 
equation for the probe soliton's amplitude has the same form as the equation for the local 
space average of energy dissipation in random cascade models in turbulence. We found 

10 



that the nth moment of the two-time correlation function and the BER contribution from 
amplitude decay exhibit power-law behavior as functions of propagation distance, where the 
exponents can be expressed in terms of the r q exponents or the Cramer function. Thus, 
our study provides a surprising and very useful perspective on the relation between disorder 
effects on weakly nonlinear systems described by perturbed NLS equations, and strongly 
nonlinear systems, such as turbulent flow. 
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